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Abstract 

The rare K + — * n + vi> decay is investigated in the context of type II two- 
Higgs-doublet model (2HDM). By using the existing experimental data of the 
branching ratio, restrictions on the free parameters of the model mji, and 
tanfi are obtained: 0.7 < tanf3 < 0.8, and 500GeF < m H < 700GeV. 







I. INTRODUCTION 



The determination of the elements of the Cabibbo-Kobayashi -Maskawa matrix 
(CKM) is still an important issue in the flavor physics. The precise determination of the 
CKM parameters will be one of the most important progresses to understand the nature, 
physics of violated symmetry. 

In this sense the rare K + — > it + vu decay has attached a special interest due its sensitivity 
for the determination of CKM parameters, in particular the element V t d, and considered 
one of the cleanest decays from a theoretical standpoint. Moreover this decay occupies 
a special place, since for this decay the short distance effects dominated over the long 
distance effects. Over the years important refinements have been added in the theoretical 
treatment of K + — > 7i + uu, long-distance contributions to the branching ratio were estimated 
quantitatively and could be shown to be essentially negligible as expected, two to three orders 
of magnitude smaller than the short distance contribution at the level of the branching 
ratio [1]. On the other hand, the calculation [2] of next-to-leading QCD correction reduced 
considerably the theoretical uncertainty due to the choice of the renormalization scales 
present in the leading order expression. Since the relevant hadronic matrix element of the 
operator S7 M (1 — 75)^7^(1 — 75)1/ can be extracted in the leading decay K + — > it°e + v. 
Conventionally, the Br(K + — > 7i + uu) is related to the experimental well-known quantity 
Br(K + — > n + e + iy) =0.0482, measured to 1% accuracy. The resulting theoretical expression 
for Br(K + — > -n + vv) is only a function of the CKM parameters, the QCD scale A^ s and the 
quark masses m t and m c . 

Experiments in the K meson system have entered new period. That the branching 
ratio of the flavor-changing neutral current (FCNC) process K + — > has been recently 

measured, and it has turned out to be Br(K + — > = (^Ijj^.lO -10 [3]. The central 
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value seems to be 4-6 times larger than the predictions of the Standard Model (SM) Br(0.6- 
1.5)xl0~ 10 [4]. Hence the rare K + — > n + uu decay is very sensitive to a new physics beyond 
the SM. Therefore, careful investigation of this decay can provide useful information about 
new physics [5]. For this reason different new physics scenarios for this decay will become 
very actual. 

In this work the decay K + — > n + uu is investigated in the framework of the two Higgs 
doublet model (2HDM). We estimate the constraints of the 2HDM parameters namely, 
tanf3 and m#, using the result coming from the measurement of [3]. Subsequently, this 
paper is organized as follows: In Section 2, the relevant effective Hamiltonian for the decay 
K + — > n + uu in 2HDM is presented. Section 3, being devoted to the numerical analysis of 
our results; and finally a brief discussion of the results is given. 



II. EFFECTIVE HAMILTONIAN 



In the Standared Model (SM) the process K + — > n + uu is described at quark level 
by the s — > duu transitions and received contributions from Z°-penguin and box diagrams. 
The effective Hamiltonian relevant to s — > duu transition is described by only one Wilson 
coefficient, and its explicit form is: 

= fas***** 9 * - ^,.(1 - 7>, (i) 

where G is the Fermi coupling constant, a is the fine structure constant, VtbV t * are products 
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of Cabibbo-Kabayashi-Maskawa matrix elements and x t — The resulting expression 
of Wilson coefficient Cn, which was derived in the context of the SM including 0(a s ) 
corrections is [6,7] 
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Here Li 2 (l - x) = J? j^dt and x^ = with /i = 0{m t ). 

At fi = m t , the QCD correction for Xi(x) term is very small (around ~ 3%), and 
rj = 0.985 is the next-to-leading order (NLO) QCD correction to the t- exchange calculated 
in [2]. 

From the theoretical point of view, the transition s — > dvv is a very clean process 
as pointed out, since it is practically free from the scale dependence, and free from any 
long distance effects. In addition, the presence of a single operator governing the inclusive 
s — > dvv transition is an appealing property. The theoretical uncertainty within the SM 
is only related to the value of the Wilson coefficient C n due to the uncertainty in the top 
quark mass. In this work, we have considered possible new physics in s — > dvv only through 
the value of that Wilson coefficient. 

In this spirit, the process s — > dvv in the context of the 2HDM has additional con- 
tributions from Z°-penguin and box diagrams through H boson exchanges. The relevant 
Feynman diagrams correspond to the transition s — > dvv has been given in [8,9]. The first 
three diagrams describe the effective Hamiltonian in the SM, while the last three diagrams 
represent the 2HDM contributions to the s — > dvv transition, due to the charged Higgs 



boson exchanges. The interaction lagrangian between the charged Higgs bosons fields and 
fermions are then given by: 

L = {2V2G F ) 1/2 \tanffU L V CKM M D D R + ctgffUnMuVcKMDL + tanffN L M E E R ] H+ + h.c. (5) 



Here, H + represents the charged physical Higgs field. Ul and Dr represent left-handed up 
and right-handed down quark fields. N E and Er are left-handed neutral and right-handed 
charged leptons. M D , Mjj, and M E are the mass matrices for the down quarks, up quarks, 
and charged leptons respectively. Vqkm is the Cabibbo-Kobayashi- Maskawa matrix, tariff- 
is the ratio of the vacuum expectation values of the two Higgs doublets in 2HDM, and it is 
a free parameter of the model. 

From eq.(5), it follows that the box diagrams contribution to the process s — > dvv 
in 2HDM are proportional to the charged lepton masses; and therefore, they are giving 
a negligible contribution. So in this model, the transition s — > dvv in eq.(l) can only 
include extra contribution due to the charged Higgs interactions. Hence, the charged Higgs 
contribution modify only the value of the Wilson coefficient Cn (see eq.(l)), and it does not 
induce any new operators (see also [8,9]): 



As we noted earlier the QCD corrections practically do not change the value of C\\. If 
so, eq. (2) and eq.(6), are plugged in eq.(l), to obtain a modified effective Hamiltonian, 
which represents s dvv decay in 2HDM: 





where x 




f- and y 




G a 
v f2 2irsin 2 6, 



w 



V td V*[X tot \s-f^(l - 7 5 )rfz/7 M (l - 75 )z/, 



(7) 



where X tot = + C?P DM . 
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However, in spite of such theoretical advantages, it would be a very difficult task to detect 
the inclusive s — > dvv decay experimentally, because the final state contains two missing 
neutrinos and many hadrons. Therefore, only the exclusive channels are expected, namely 
K + — > 7T + z/z/, are well suited to search for and constrain for possible "new physics" effects. 

In order to compute K + — > 7i + uu decay, we need the matrix elements of the effective 
Hamiltonian eq.(7) between the final and initial meson states. This problem is related 
to the non-perturbative sector of QCD and can be solved only by using non-perturbative 
methods. The matrix element < n + | H e ff | K + > has been investigated in a framework 
of different approaches, such as chiral perturbation theory [10], three point QCD sum rules 
[11], relativistic quark model by the light front formalism [12], effective heavy quark theory 
[13], and light cone QCD sum rules [14,15]. As a result, the hadronic matrix element for the 
K + — > -k + vu can be parameterized in terms of form factors: 

< 7i | s 7/1 (l - lb )d | K > = fl + {q 2 ){pK+Pn)^ + f-q„, (8) 



where = px — Pw, is the momentum transfer. In our calculations the form factor /_ 
part do not give any contributions since its contribution ~ m u =0. After performing the 
mathematics and taking into account the number of light neutrinos N v = 3 the differential 
decay width is expressed as: 

Jq " 2%W^ + 1 VtbVts 1 A (1 ' r+ ' S) 1 Xtot 1 1 /+ {q ' 1 ' (9) 



m? , , „2 



where r+ = -^r- and s = -% — . 

Similar calculations for K + — > -n°e°u lead to the following result: 

" r(A 'V/ e+P> — IK., I 2 A^ (1 , r ., s) < + 1 |», (io) 
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where r_ = — f 2 -, and A(l,r ± , s) = 1 + r\ + s 2 — 2r±s — 2r± — 2s is the usual triangle 

m K+ 

function. In derivation of eq.(10) we neglect the electron mass, and we use the form factors 
f+ + = \f2ff which follows from isotopic symmetry. Using eq.(9) and eq.(10) one can 
relate the branching ratio of K + — > 7i + uu to the well known measured decay K + — > n°e°i> 
branching ratio: 

B(B + -> K+vv) = k 

where k = r K+ ** B & + r?? e+v) A 8 = 4.11. 10" 11 . 

Here r K + = 0.901 summaries the isospin-breaking corrections which come from phase 
space factors due to the difference of masses of n + and n°. 

In derivation eq.(ll) we have used the wolfenstein parametrization of the CKM matrix, 
in which each element is expanded as a power series in the small parameter 
A =| V us |=0.22, Aj = V* s V i( i and P(K + ) represent the sum of charm contributions to the 
two diagrams including the (NLO) QCD corrections [2]. At m c = 1.3 GeV, A^ s = 0.325 
GeV and at renormalization scale \i c = m c in [16] it is found that Po(K + ) = 0.4 ± 0.06. 

III. NUMERICAL ANALYSIS 

In the numerical analysis, the following values have been used as input parameters: 
G F = 1.17 .10~ 5 GeV~ 2 , a = 1/137, and A = 0.22 As we noted early we used the wolfeustein 
parametrization of CKM matrix elements. In this parametrization 
ImXt = A 2 X 5 r], Re\ c = -A(l - A 2 /2), and ReX t = -A 2 X 5 (1 - p). 

The parameter A determines from b — > c transition and its A = 0.80 ± 0.075 [17]. The 
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other two CKM parameters p and rj are constrained by the measurements of | V u b/V c b |, 
x<i(B® — B®) mixing, and | e | (the CP violation parameter in the kaon system). For typical 
values of the necessary input parameters of p and i] we have adopt the following two sets : 



set I 



p = 0.06 
r] = 0.35 



set II 



p = -0.25 
77 = 0.3 



(12) 



The free parameters of the 2HDM model which we have used namely tan/3 and m H are 
not arbitrary, but there are some constraints on them by using the existing experimental 
data. These constraints are usually obtained from B° — B°, K° — K° mixings, b — > 57 decay 
width, Rb = Y{z^h2drlns) ' an d semileptonic b — > cu T r decay which are given by [18] as 

0.7 < tan/3 <0.6( I ^ ? ), (13) 

where as a lower bound for the charged Higgs mass m H > 300 GeV at p = 5 scale has been 
estimated in 2HDM [19]. If these constraints are respected, an upper and lower bound for 
ctgf3 is extracted: 

0.004 < ctgp w 2. (14) 



In Figures 1 and 2, we represent the branching ratio of K + — > ix + vh> as a function of ctg(3 
for various values of m#, and as a function of m# for various values of ctgf3. For illustrative 
purposes we consider three values of ctg/3, namely ctgf3 =1, 1.5 and 2 and we allow m H to 
range between 300 GeV and 1000 GeV, and then we consider three values of uih, namely 
m#=300, 500, 1000 GeV and we allow ctg/3 to range between to 2. It can be seen that 
for ctgf3 = 1, the branching ratio (BR) for K + — > -n^vv decay increases slowly with the 
increasing of m H ; whereas, for larger values of ctg/3, the BR decreases at all values of m#. 
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Furthermore, when the ctg(5 is increased the BR rapidly grows up. Therefore, it can be 
concluded that the main contribution to the decay width comes from the charged Higgs 
exchange diagrams (see [8,9]). 

The question now is; what kind of restrictions on tanf3 and m# can be obtained if the 
recent experimental result of Br(K + — > ix + h>h>) = (4.2l|; 5 ). 1CT 10 [3] is respected that is: 

(0.7 < BR exp - < 13.9). 1(T 10 , (15) 

and whether or not it coincide with the restrictions given in [18]. For this aim, in Figure 3 
we present the dependence of tanf3 on m# using both sets of values of p and rj. We see that 
when 300GeV < m H < ITeV it gives: 

0.18 < tanp < 0.5 ±0.2 (setl), (16) 
0.18 < tan(3 < 0.8 ±0.3 (setll). (17) 

If we use the lowest bound for tanf3 = 0.7 (see eq.(13)) we see that the set I predictions is 
ruled out and for set II we have small room for tan(3, namely from eq.(16) and from eq.(17) 
we have: 

0.7 < tan(3 < 0.8. (18) 

If we increase a little bit the upper bound and if we put a lower value for m#=500 GeV we 
can see that in this case 

0.7 < tan(3 < 0.9. (19) 



Using these results we can conclude that the mass of the charged Higgs boson must be 
lie in the interval 
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500GeV <m H < 700GeV. 



(20) 



In conclusion, using the experimental result of the branching ratio for K + — > tt + vv and 
the CLEO measurements on b — > 57 [20] we find new restrictions on the free parameters 
tanf3 and ran of the 2HDM model. In summay it is found that the contribution of type 
II two-Higgs-doublet model to the branching ratio is exceed at most by ~ 20% from the 
standard model ones. 
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Figure Captions 



Figure 1 : The dependence of the Br(K + — > ix+vv) on ctgfi at fixed values of m H . 
Figure 2 : The dependence of the Br(K + — > -k + vu) on m H at fixed values of ctgfl. 
Figure 3 : The dependence of tan/3 on the charged Higgs boson mass m H - Curves (A, B), and 
(C, D) describes upper and lower bound of the experimental values of the Bi(K + — > -n + vv) 
for Set I and Set II values of p and n respectively. 
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